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Abstract 



Numerical simulations, based on the core-nucleated accretion model, are presented for 

the formation of Jupiter at 5.2 AU in 3 primordial disks with three different assumed 

p^ ' values of the surface density of solid particles. The grain opacities in the envelope of the 

protoplanet are computed using a detailed model that includes settling and coagulation 

of grains and that incorporates a recalculation of the grain size distribution at each point 

in time and space. We generally find lower opacities than the 2% of interstellar values 

O ' used in previous calculations [Hubickyj, O., Bodenheimer, P., Lissauer, J. J., 2005. Icarus 

fi ! 179, 415-431; Lissauer, J. J., Hubickyj, O., DAngelo, G., Bodenheimer, P., 2009. Icarus 

^ ' 199, 338-350]. These lower opacities result in more rapid heat loss from and more rapid 

contraction of the protoplanetary envelope. For a given surface density of solids, the new 

calculations result in a substantial speedup in formation time as compared with those 

previous calculations. Formation times are calculated to be 1.0, 1.9, and 4.0 Myr, and 

1^ ■ solid core masses are found to be 16.8, 8.9, and 4.7 M^, for solid surface densities, a, 

f~^ , of 10, 6, and 4 g cm~^, respectively. For a = 10 and cr = 6 g cm~^, respectively, these 

OO ■ formation times are reduced by more than 50% and more than 80% compared with those 

in a previously published calculation with the old approximation to the opacity. 



f~^ ' Keywords: Planetary formation, Jovian planets, Jupiter: interior, Accretion 
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1. Introduction 



^ ' The c ore-nucleated ac cretion theory of giant planet formation was effectively origi- 

5t , nated bv ISafronovl (J1972I) . who realized that a major part of the process involved the 



accumulation of solid bodies. Consider the accretion of relatively small objects (plan- 
etesimals) onto a significantly larger embryo. If Mcore is the mass of the embryo, then 
the fundamental equation for its growth, in the absence of gas, is 

^^ = .Rl^^^a^F, « .i?Lc^f^[l + i^f] , (1) 
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where ttR^^^^ is the effective geometrical capture cross-section, a is the surface density 
of soUd material (planetesimals) , v^ is the escape velocity from the embryo, Q, is the or- 
bital frequency, and v is the relative velocity of embryo and planetesimal. The quantity 
in brackets is t he gravita tional enhancement factor Fg, here given in the (2 -|- 2)-body 
approxim a tion. ISafronovl mentioned the later capture of gas to form giant planets, and 
CameronI (|l973r) stated that Jupi ter could form by gas acc retion once a solid core of 



about 10 M0 had been accreted. IPerri and CameronI ( 19741) then constructed an adia- 
batic, hydrostatic equilibrium model of a protoplanet with a core and a gaseous envelope, 
extending outward to the Hill radius. They were able to show that there was a critical 
mass for the core, above which the envelope was unstable, implying that rapid gas ac- 
cretion could occur. However, the value of this critical core mass was at least 70 M^, 
much higher than the deduced core m asses of Jupiter a nd Saturn. The assumption of 
adiabaticity turns out to be incorrect (Steve nsonl Il982 ). as more detailed calculations 
show that radiative transport with subadiabatic gradients is able to carry the energy 
ou twards in signifi cant portions of a growing giant planet. 

Mizund ( 19801 ) constructed models of a protoplanet including both radiative and 



convective energy transport in the envelope, along with an energy source provided by 
the planetesimals accreting onto the core at a constant rate, typically Afcorc = 10"® 
Mg/yr. The model assumed strict hydrostatic and thermal equilibrium, and the outer 
radius was assumed to be the protoplanet's Hill radius. Again a critical core mass was 
found, above which the envelope was unstable, and the value was about 12 M0, close 
to the estimated core masses of the giant planets. Furthermore, the critical value was 
practically independent of the distance to the Sun, partly as a result of the assumption 
of the same constant value of Afcoro at all distances. 

Th e assumption of strict hydrostatic equilibrium was relaxed bv lBodenheimer and Pollack 
( 19861 ). who allowed quasi-static contraction and calculated the first evolutionary se- 
quences applied to the core-nucleated accretion model. The planetesimal accretion rate 
was assumed to be constant, also at 10~® M^/yr. They found that gravitational con- 
traction (but not dynamical collapse) began to become important at about the point 
where core and envelope masses were equal, at a value which later came to be known as 
the crossover mass. Rapid gas accretion also began near the crossover mass, which they 
found to be about 16.5 M®. 

The next major step was a series of calculations bv lPollack et al.l ( 19961 ) who allowed 

a variable core accretion rate given by eq. ([T]). The simulations contained thre e major ele- 

ments: (1) three-body accretion cross-sections of solids onto the embryo ( Greenzweig and Lissauer . 



19921 ) to determine F^, (2) a modified stellar evolution code (jHenvev et al.lll964l ) to follow 



the structure and evolution of the gaseous envelope, and (3) a cal culation of the trajecto- 
ries of planetesimals through the gaseous envelope (jPodolak et al.. 1988 ). to deduce the 
deposition of solid mass and energy in the envelope and to determine the protoplanet's 
effective capture radius, taking into account gas drag. The often-quoted result of these 
calculations is that it took Jupiter 8 Myr to form in a solar nebula with cr = 10 g cm~^, 
three times the density in the minimum mass solar nebula, under the assumption that 
the opacity due to small grains in the cool outer envelope was provided by grains with 
an interstellar size distribution and with interstellar abundances with respect to the gas. 
It was also shown that if the grain opacities were arbitrarily reduced to 2% of interstellar 
values, the formation time was reduced to about 3 Myr. 

Subsequent observations of protoplanetary disks and modelling of the interior of 
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Jupiter revealed two major problems with the iPollack et alJ ()1996ll baseline-case results 
for the formation of Jupiter. First, that case gave a formation time longer than the mean 
lifetime of protostellar disks, and, second, it produced a core mass of about 16 M^, above 
the plausible range at the time of (0-1 1) Mm, derived from matching interior models of 
the present Jupiter wi th observations (jSaumon and Guilloti . 120041 ) . Further calculations 



(IHubickvi et all 120051 ) considered, with generally improved physics, the effects of the fol- 



lowing parameters: (1) the opacity arising from small grains in the envelope, (2) the solid 
surface density a in the protoplanetary disk, and (3) a possible cutoff in core accretion 
rate as a result of competition for solid particles by ne ighboring embryos. The reduction 
in opacity was justified by calculations ( Podolakll2003l ) showing that grains entering into 



the protoplanetary envelope would rapidly coagulate, settle, and eventually evaporate, 
resulting in an actual opacity well below interstellar values. In the actual protoplanetary 
formation calculation the grain opacity was simply reduced to 2% of interstellar values, 
with the result that the formation time was reduced to under 3 Myr, half that found for 
interstellar values of grain opacity. The final core mass, however, still remained near 16 
M^. A reduction in a from 10 to 6 g cm~^, again with the opacity at 2% of interstellar, 
increased the formation time to 13 Myr and reduced the core mass to 8 M^. It turned 
out that the constraints of short formation time and small core mass could be satisfied 
simultaneously only if the accretion rate of solids into the protoplanet were arbitrarily 
cut off at some time, simulating roughly the effect of planetesimal accretion by neighbor- 
ing protoplanets. A cutoff at Mcorc — 5 M^ led to a formation time of 4.5 Myr, while a 
cutoff at 10 Mgg gave a formation time o f 1 Myr. 

In the planet formation simulations of lPollack et al.l ( 19961) and lHubickvi et al.l ( 20051) . 



simplified surface boundary conditions were applied during the rapid gas accretion phase. 
First, the gas density in the protoplanetary disk was assumed to be constant in time. 
Second, the gas accretion rate was capped at 10~^ M^/yr, a rough estimate of the rate 
at which the disk could supply gas to the protoplanet as a result of its viscous evolution. 
Third, the out er boundary of the p rotoplanet was set at the modified accretion radius 
Ra, defined bv IBodenheimer et al.l (|200Q ') to be 



Cs/'«l-t- (k2RH) 

where Rh is the Hill sphere radius, Cg is the sound speed in the disk, Mp is the planet's 
mass, and fci, fc2 are constants, both set to unity. Equation ^ is an interpolation 
between the Bondi accretion radius and the Hill sphere radius, such that when Rh is 
large, Ra reduces to the Bondi radius, and when Rh is small, Ra reduces to the Hill 
radius. 

These assumptions were modified bv iLissauer et al.l (|2009() to take into account the 



results of new three-dimensional numerical simulations of protoplanetar y disks accreting 
onto p rotoplanets. These simulations employed the code developed by iD'Angelo et al.l 



(|2003f ). The results from the 3-D simulations showed that ^2 ~ 0.25, that is, the planet 



can collect material from the disk only within a quarter of its Hill radius, at best. The gas 
accretion rates from the 3-D simulations, which included the effects of gap opening, were 
used whenever they were less than the accretion rate needed to maintain the boundary 
of the planet at Ra', thus the arbitrary cap on this rate was removed. The parameters 
investigated in this set of simulations included (1) the size of the region from which the 
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protoplanet can accrete, (2) the viscosity in the protoplanetary disk, which has a strong 
influence on the disk accretion rate onto the protoplanet once the envelope has shrunk 
and the planet is massive enough to begin clearing gas from nearby its orbit, and (3) 
the time-dependence of the gas density in the protoplanetary disk in the vicinity of the 
protoplanet. The grain opacity was set to 2% interstellar and cr = 10 g cm^'^. The 
results of these improved simulations showed that Jupiter's formation time was still less 
than 3 Myr, and that a relatively low-viscosity disk (a ~ 4 x lO""*) is consistent with 
the formation of planets of about 1 Jovian mass, while higher viscosity produced planets 
that were more massive than most observed extrasolar giant planets. 

In all of the s imulations just men t ioned, as well as in detailed formation simulations 
by other groups (JAlibert et all . l2005t iFortier et al.l . l2007t llkoma et all . 120001 ) . the opac- 



ities arising from grains in the protostellar e nvelope have been approximated either by 



interstellar grain values ([Pollack et al.l 119851 ) or by a reduction of those values b y some 



arbi trary amount, typically a facto r 50. However, the simulations of IPodolakI (2003) 



and iMovshovitz and PodolakI ( 20081 ) . in which grain coagulation and sedimentation were 



calculated for a few particular protoplanetary envelopes, showed that the ratio of actual 
grain opacity to interstellar grain opacity was highly variable, depending on the depth in 
the envelope. Near the surface of the planet the opacities were close to interstellar, while 
deep in the envelope they were reduced by a factor of up to 1000. Those opacity calcula- 
tions, however, were not coupled to the evolution of the protoplanet, and they included 
only the outer radiative layer in the atmosphere, which typically has a temperature of 
600 K at its inner boundary. 

In the present paper we ex tend the basic grain settling and coagulation model of 
Movshovitz and PodolakI ( 20081) by including the entire depth of the envelope down to a 



temperature of 1800 K, where practically all grains have evaporated. Since some layers of 
the outer part of the envelope are convective, we also modify the grain simulation to take 
into account convective effects. We then couple the grain calculation to the evolution 
of the protoplanet, recalculating the grain size distribution and opacity in every layer 
of the protoplanet at every time step of the evolutionary calculation. These significant 
improvements in the core accretion model allow us to investigate the following questions: 
(1) How is the formation time for a Jupiter-mass planet at 5.2 AU modified compared 
with the earlier approximate calculations regarding the opacity? (2) Most of our previous 
simulations for the formation of Jupiter required an initial solid surface density of 10 g 
cm~^ in the protoplanetary disk. How far can this assumed density be reduced such 
that Jupiter still forms within a dissipation time of the disk - roughly 3 Myr ~ and 
what is the resulting solid core mass? The following sections describe the nature of the 
grain simulations, the basic idea of the core-nucleated accretion model, how the two 
calculations are integrated, and our new results for four formation runs, up to the onset 
of rapid gas accretion, for a Jupiter-mass planet at 5.2 AU. 

2. Grain Physics 

This section describes, for a protoplanetary envelope with a given distribution of 
temperature, density, and convective velocity as a function of distance to the center, how 
grain settling and coagulation modify the grain size distribution and opacity at each level. 
The size distribution of the grains as a function of height is determined as follows: At 



each atmospheric level the number of grains per unit volume, n{m, r, t), evolves according 
to 



dn{m, r, t) 
dt 



Q{m,r,t) 



dn{m,r,t) 
dt 



+ 



dn(m, r, t) 



dt 



(3) 



where the derivatives are Eulerian. Here m is the mass of the grain, r is the distance 
from the center of the planet, and t is the time. 

The first term on the RHS of Eq. ([3]) is a source term, and refers to those grains that 
are carried into the envelope either by the accreted gas or by planetesimal ablation. We 
do not consider the effects of evaporation and condensation on the grain size distribution. 
A grain is instantaneously evaporated when its temperature reaches f 800 K. (We assume 
that the grains are composed of silicates. The error introduced by this assumption is 
probably small compared with the overall uncertainty in the physics of grain structure 
and agglomeration. This point is furthered discussed below.) 

The second term on the RHS of Eq. ^ refers to the evolution of the grain size 
distribution by coagulation. We treat this by solving the Smoluchowski equation 



dn{m, r, t) 
di 



1 



K{m' , m — m' , r, t)n{m' , r, t)n{m — m', r, t)dm' 



K{m, m , r, t)n{m, r, t)n{m , r, t)dm . 



(4) 



The coagulation kernel, K{m,m' ,r,t) is essentially the probability of a collision be- 
tween grains of mass m and m' that results in the two grains sticking t ogeth er to 
form a new grain. The details of how K is computed are given in iPodolakI (|2003f ) and 
Movshovitz and Podola^ (|2008[ ). 

In practice, the grain size distribution is approximated by a finite number of size bins, 
where the zth bin represents a grain of radius 



2^/3 



flo 



and g p is a scaling constant, set in this case to 1 jira. It has been shown (jMovshovitz and Podolak , 
2008! ) that if the minimum grain size in a protoplanctary envelope is reduced to 0.1 //m, 
the results are practically identical to those for a minimum size of 1 /im, because of rapid 
coagulation of the small particles. The number of bins is chosen so that the size range of 
interest can be adequately represented. As a result, the integrals in Eq. ^ are replaced 
by sums. 

The third term on the RHS of Eq. ([3]) is a transport term. In the absence of convection, 
the transport is given by 



dn{m, r, t) 



dt 



-V- [n{m, r, t)vsed{m, r, t)] 



(5) 



Here Vsed is the sedimentation velocity of the grain, and it is a function of the grain size 
and t he local values of gravity, gas density, and temperature (jMovshovitz and Podolakl 
2008h . 

If convection is present, it affects the results in two ways. First, small grains that 
are entrapped by a convective eddy reach higher speeds than normal, while large grains 
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are less affected by these eddies. This effect causes the relative v elocities between smah 
and large grains, and the resultant coagulation kernels, to change. IWeidenschillina (|1984l ) 
gives expressions for grain relative velocities in the limits where both grains are large, one 
is large and one is small, and bo th small (his Eqs. 15, 16, and 18, respectively). More 
recently, lOrmel and Cuzzil (J2007l) give an expression for t he relative veloc ity between 
particles of arbitrary size (their Eq. 16), that matches the IWeidenschillinsJ expressions 
in the relevant limits. We use the lOrmel and Cuzzil expression, modified to include 
the effects of the systematic (sedimentation) velocity of the grains. The sedimentation 
velocity shifts the boundary between small (coupled to the gas) and large (uncoupled) 
grains, and the difference in sedimentation velocity between two grains is also added in 
quadrature to the relative velocity due to turbulence. 

In practice, however, the enhancement of collision probability due to turbulence, 
although fully included in the calculations, is not terribly important. For the most part, 
relative velocities between grains of different sizes are determined by the difference in 
sedimentation speed, and the relative velocities between small particles of similar size 
are determined only by their random thermal motion, because they are coupled to the 
gas. The turbulence does add collisions between large grains of equal size. These grains 
are too large for significant thermal motions, and with no difference in sedimentation 
speed they have no appreciable relative velocities, except those induced by turbulence. 
Because this enhancement is only apparent for large grains (« 3 x 10"^ — 1 cm), the 
effect on the size distribution and opacity is negligible. 

A second e ffect of convection is that the transport term becomes more complicated. 



PodolakI (|2OO30 treated convective transport with an eddy diffusion approximation. This 
leads to numerical difficulties and requires very short time steps. As a result, it is un- 
suitable for use in the current evolutionary models. We have therefore treated convective 
transport with a simplified algorithm that retains the major characteristics of convective 
mixing, explicitly conserves mass, and is very computationally efficient. For each layer 
in the envelope and for each grain size, we decide if convection is important or not based 
on whether two conditions are satisfied. The first is that the time it takes for a grain to 
become entrapped in the largest eddy be shorter than the lifetime of that eddy. If the 
grain is small enough so that the time to entrap it is short, the grain responds to the 
convective motion. For this calculation, we have taken the length of the largest eddy to 
be 1.5 times the pressure scale height. The eddy lifetime was taken to be this length 
divided by the speed of the eddy. The second condition that is relevant to convective 
transport is that the convective speed of the eddy be greater than the sedimentation 
speed of the grain. If this condition is not met, the grain sediments faster than the 
convection can mix it. Going back to the first condition, note that there are actually two 
requirements for a grain to be carried by an eddy. First, as stated above, the stopping 
time must be shorter than the eddy turnover time. Second, the stopping time must be 
shorter than the grain crossing time (eddy size/grain velocity). In our case, the crossing 
time is ~ L/vsed, and the turnover time is ^ L/vconv, where L is the eddy size. Since 
the second condition requires v^onv > Vsed, the crossing time is longer than the turnover 
time, which wc have already required to be longer than the stopping time. Thus the 
requirement (stopping time < crossing time) is automatically met. 

For each layer in the envelope and for each bin size we determine whether or not these 
two criteria are met. If they are, convection will be the dominant mode of transport. 
We then determine, for each grain size, which regions (i.e., sets of contiguous layers) of 
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the envelope will mix that grain size. In those regions we determine the total number of 
grains, and redistribute them so that the ratio of mass density of grains to mass density 
of gas is constant. Grain sedimentation velocities range from millimeters to meters per 
second. We find that almost all grains are mixed for convective velocities larger than 
~ 1000 cm s-^ 

This algorithm has the disadvantage that it allows for either complete convective 
mixing or no convective mixing. Partial convective mixing is not considered. However, 
it has the advantage that it explicitly conserves the grain mass in the region, and it 
allows much larger time steps than the eddy diffusion algorithm. Numerical tests show 
that the resultant grain size distributions behave in accordance with what is expected for 
convective mixing, so this seems to be a useful approximation to convective transport. 

Once the grain size distribution is established, the opacity is compute d using an ap- 



proxi mation to Mie scattering suggested by Cuzzi (personal communication jMovshovitz and Podolak . 



20081 ). The effective cross-sections for absorption and scattering are calculated for each 
grain size and frequency, and combined to produce the opacity K^{ai). This function 
is then integrated over the grain size distribution and over the blackbody spectrum to 
produce the Rosseland mean opacity for each level in the atmosphere. 

We have assumed that the flux of grains into the uppermost layer of the envelope is 1% 
of the gas accretion mass flux, given a near-solar composition in the protoplanetary disk. 
This will cause us to overestimate the opacity for the following reason. It assumes that 
all of the condensed material in this region is in the form of small grains. This is certainly 
not true, since many of the grains will have accreted into larger planetesimals. These 
planctesimals are eventually accumulated by the planet, but their grains are released 
into the envelope at much deeper levels and thus have much less of an influence on the 
opacity of the outer layers. 

Our model also assumes that all of the grains are composed of silicates. In fact, 
they are likely to be composed of a combination of silicates, organic material and water, 
in roughly equal fractions by mass. This has two major consequences. First, although 
the silicates and organics will remain solid at the temperatures relevant to the outer 
envelope, the water will evaporate. This means that we have overestimated the mass of 
grain material at all levels hotter than ^ (170 — 200) K and the true grain opacity will 
be correspondingly lower. 

A second effect is that the bulk density of the grains in these upper layers will be 
lower than we have assumed, and this means that the grains will sediment more slowly 
and their number density will be higher than we have computed. As a result there should 
be some increase in the opacity. Below we estimate the magnitude of this effect. 

In the uppermost layers of the envelo pe, there is very li ttle grain growth and most 



of the grains are the same size (jMovshovitz and PodolakLl2008i) . To first approximation. 



the time for such a grain to settle out of a layer of thickness I is 

tsed = l/Vsed (6) 

where Vsed is the sedimentation speed. The time between grain collisions is roughly 

t -- (7) 

^coaq — ■ \* } 

Vd 

Here A is the mean free path between grains, and Vd is the thermal speed of a grain. 
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^itcoag is short compared to tsed, the grains will collide and grow before they sediment 

trnnn- The mean free 



out of the region. The grain size can be estimated by setting tg^d 
path between grain collisions will be 



coag 



A = 



1 



V2(2a) 



(8) 



where n is the number density of grains, and a is the grain radius. The Knudsen number 
in this region is large so that the sedimentation speed can be written as 



Vsed 



PgVg 



-ga, 



(9) 



where g is the acceleration of gravity, pd is the bulk density of a single grain, pg is the 
mass density of the gas, and Vg is the thermal speed of a gas molecule. The mass flux 
through the layer is denoted by F, where 



An ., 

F = —npdVseda ■ 

The thermal speed of a molecule (or grain) of mass m is given by: 



(10) 



8A;T 



(11) 



where T is the temperature and k is Boltzmann's constant, 
using the above relations, we find 



3VT2kfFlp 



>l 



7/2 2 

T'Pd r 



2/9 



Setting tsed = tcoag and 



(12) 



The cross-section for scattering or absorption is the geometric cross-section, tto^, multi- 
plied by some eflficiency factor, Q, that must be computed from Mie theory. Typically, 
Q « 1 for grains much smaller than the wavelength of the impinging radiation, and 
roughly equal to one for larger grains. Since we are interested in Rosseland mean opac- 
ities at temperatures of around 100 K, the wavelengths will be in the 20 — 30 micron 
range and the grains are in the 1 micron range in the upper levels of the envelope. The 
corresponding size parameter is 

27ra 
X-—, (13) 

where A is the wavelength of the radiation. 

For s mall x the extinction is mostly due to absorption, and a useful approximation 
for Q is (JMovshovitz and Podolafl l2008l ): 



Q 



2Axnrni 



("r 



1 



(14) 



where rir and Ui are the real and imaginary indexes of refraction, respectively. Typically 
rir ~ 1.5 and rii w 0.1 so Q ~ 0.2a;. In our case, x « 0.25, so Q « 0.05. 



The opacity per gram of gas is 

Qna'^n 3QFvg 



Ap^^ga^ 



(15) 



All of these quantities are interconnected, since the grain size a is also a function of grain 
density to the power —(7/9) . As a result, the opacity in these upper layers is expected 
to depend only weakly on the assumed grain density. Thus, even if the grains are made 
of pure silicates with a density of ~ 3 g cm^^ or of pure ice with a density 1/3 as large, 
the opacity should change by roughly a factor of 1.6. In practice the change should be 
even smaller. 

3. Core Accretion Model 

The physics of t h e core-nucleated accretion model are for the most part described by 
Pollack et al.l (|l996l) , iBodcnhcim er et al.l (|2000l ) , iHubickvi et al.l (|2005l ) , and lLissauer et al. 



( 20091 ) . The calculation employs the three basic elements mentioned in Section [TJ The 



first element is essentially Eq. ([T|), which is used to cal culate the accretion rate o f solids . 
The value of Fg is obtained from the numerical fits of iGreenzweig and Lissaueii ( 1992I) . 



a 

assuming a planetes imal radius of 1 km . Other more detailed calculations of the 
cor e accretion rate (IKokubo and Idal . 120001 ) show in general that the rates calculated 



by IGreenzweig and Lissaueij (|1992 i) should be reduced because they do not completely 



take into account the stirring and increase in relative velocity v of the planetesimals as 
a result of gravitational interactions with the embryo. However, the accretion rate for 
smaller planetesimals is faster than that for 100 km. Thus our accretion rate should be a 
reasonable approximation to the actual rate assuming that the planetesimals determining 
the velocity distribution are somewhat smaller than 100 km. 

As in previous calculations, we take the feeding zone from which the embryo can 
accrete planetesimals to extend 4 Hill sphere radii on either side of the orbit, and assume 
that the surface density a is spatially uniform within that zone. As the planet accretes 
material, a in the feeding zone is recalculated, again assuming uniform density, taking 
into account depletion of planetesimals by accretion onto the embryo and expansion of 
the feeding zone into undepleted regions as the embryo's mass increases. 

The second element of the code calculates the interactions b etween planetesimals and 
the gaseous envelope of the protoplanet as they fall through it (Podolak et al., 1988). The 
equations of motion of the planetesimals are integrated, for a variety of assumed impact 
parameters, taking into account the gravitational fields of the core and the envelope 
as well as gas drag. The critical impact parameter inside of which the planetesimal is 
captured is determined, and the effective capture radius i?coro (Eq. H]) is assumed to be 
the periapsis altitude of the critical orbit. This radius can be several times larger than 
the actual solid core radius, once a significant gaseous envelope has accreted, so this 
effect is important in speeding up the accretion rate of solids. 

For those planetesimals with impact parameters less than the critical value, trajec- 
tories are calculated to determine the deposition of mass and energy into the envelope. 
Effects included are ablation and evaporation of planetesimal material, heat associated 
with phase changes, and fragmentation when the dynamical pressure of the gas exceeds 
the compressional strength of the planetesimal. During the earliest phases, when the 
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envelope has low mass, planetesimals plunge all the way to the core. At later stages 
they can be fully ablated or fragmented in the envelope. It is assumed that the remains 
of the planete s imal later sink to the core, releasing gravitation al energy in the process 



( Pollack et al.Ul996l ). This assumption is not entirely accurate: llaroslavitz and Podolak 



20071 ) have made a calculation of how much of the heavy-element material should ac- 



tually dissolve in the envelope. Ice should be almost entirely dissolved, but most of the 
organics and rock would sink to the core or to the layers just outside it. They conclude 
that the overall evolution would not be substantially changed if this effect were included; 
the accretion of the protoplanet would be slightly sped up and the core mass would be 
reduced to about 2/3 the calculated value, with little change in the overall metal content. 

The third element of the simulation is the solution of the four differential equations 
of stellar structure for the gaseous envelope, with the nuclear energy term replaced by 
the energy delivered by accreting planetesimals and with the gravitational energy terms 
included. The adiabatic temperature gradient is assumed in convection zones. At the 
core-envelope interface the radius is set to that of the outer edge of the core, calculated 
from Mcoro and the assumed core mean density of 3.2 g cm~^. Also, the luminosity 
Lr is set to the energy deposition rate for the planetesimals that hit the core. Outside 
the core, the energy supplied by ablated and fragmented planetesimals is included as a 
source term in the energy equation, along with the energy derived from gas accretion 
and compression. 

At the surface, mass is added at a sufficient rate to mai ntain an outer radius Rjy = 
Ra, and Eq. ^ is used with fc2 — 0.25, as reported in iLissauer et al.l ( 20091 ) from 



the results of three-dimensional calculations of disk-planet interaction. Otherwise the 
density and temperature at the surface are set to assumed nebular values Pnob, T,^ch, 
respectively. These values are constant in time. The simulations are stopped when the gas 
accretion rate reaches the disk flow limit dMn^/dt, as de termined from three- dimensional 
calculations with disk viscosity parameter a — ix 10~^ ( Lissauer et al.l . 120091 ). The lower 
curve in figure 3 of that paper gives the limiting gas accretion rate as a function of planet 
mass. In the various cases of the present calculation, dMn^/dt w (1.5 to 3) x 10^^ M^ 
yr^^ when the limit is first reached. 

The equation of state of the gas is taken from ISaumon et al.l ( 1995[ ). interpolated 
to our assumed composition of hydrogen mass fraction X = 0.74, helium mass fraction 
Y = 0.243, and metal mass fraction Z = 1 — X — Y = 0.017. Although the equation of 
state in the outer, low-density layers is essentially that of an ideal gas, the inner regions 
near the core can be significantly non-ideal once the envelope has acquired sufficient 
mass. 

The Rosseland mean opacity calculation has three components. At temperatures 
above 3500 K the molecular opacities of [Alexander and FergusonI ( 1994 ) are used. In 
practice the details of the opacities in this region are unimportant because the energy 
transport is almost always by convecti on. In the te mperature range 200-3500 K the 
molecular opacities, without grains, of Freedman et al.l ( 2008) are used. The grain opacity 
calculation described in scction[2]is then included in the temperature range 100-1800 K. 
In the higher end of this range the grain opacity can be so low at times that the molecular 
opacity dominates. Grains are assumed to be completely evaporated above 1800 K. 

The grain calculation and the stellar envelope calculation are combined as follows. 
After each evolutionary time step At for the envelope, grains are added at the surface in 
accordance with dMonv/dt, assuming a dust-to-gas ratio of 0.01. Grains are also added 
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where the ablation/fragmentation calculation for the incoming planetesimals deposits 
them, typically fairly deep in the envelope. For the grains, 34 size bins are used, with the 
smallest at 1.26 /zm and the largest at 2.58 mm. Incoming grains, represented by the term 
Q(m, r, t) in Eq. ^, are always placed in the smallest bin. The envelope calculation then 
provides to the grain calculation a table of temperature, density, and convective velocity 
as a function of radius, as well as the core and envelope masses. Velocities are estimated 
by the stellar structure code through the mixing-length theory. An arbitrary number of 
distinct convection zones is allowed. 

The grain evolution, involving settling and coagulation, then proceeds in the region 
from the surface down to 1800 K. A large number of sub-timesteps is used, adding up to a 
total time of At, as the number of grains in any size bin is not allowed to change by more 
than 4% per sub-timestep. The number of sub-timesteps can be 10,000 to 100,000, so 
this part of the calculation is quite time-consuming. Once the grain calculation has been 
completed, with a new opacity provided at each depth, the opacities are fed back into 
the envelope calculation in the form of a table of Rosseland mean opacity as a function 
of mass fraction, integrated inward from the surface. Mass fraction is used rather than 
temperature or density because, first, the temperature is practically isothermal in the 
outer layers, and, second, temperature and density values can fluctuate somewhat as a 
result of the numerical procedure for adding mass to the envelope. The current grain 
size distribution is stored at each depth for use in the next time step. The envelope 
calculation then proceeds to the next model. Once the model has been completed, the 
Lagrangian zones have moved in radius with respect to the previous model. The grain 
number densities, tabulated as a function of radius, are interpolated onto the new grid 
by a procedure that conserves the total number and mass of grains. The results for 
the grain opa city are qualitatively quite s imilar to those obtained for specific envelope 



models bv Movshovitz and PodolakI (|2008l ). In the outer zones, where small grains are 
continuously added from the protoplanetary disk, the opacity generally falls in the range 
0.1-1 cm^ g~^j not too different from opacities derived from grains with interstellar 
size distribution and abundance. In contrast, at the deeper levels the typical grain size 
increases sufficiently so that the mean opacity can drop by up to 4 orders of magnitude 
relative to interstellar values! 

4. Calculations and Results 

The planet initially consists of a core of 1 M^ and an envelope of about 10~^ M^. 
The protoplanet is located at 5.2 AU in a solar nebula disk, with the solid surface density 
a treat ed as a parame t er. Based on our pas t calculations for the time required to reach 
1 Me (JHubickvi et all . l2005t iLissauer et all . [2009), the starting time is set to 2.8 x 10^ 



yr, 4.85 x 10^ yr, and 7.28 x 10^ yr for a = 10, 6, and 4 g cm~^, respectively. No orbital 
migration is included. The quantity Tnob is set to 115 K, and p^oh = o'xv/i'^H), where 
axY = 70(7. The scale height of the gas in the disk H = 0.05a, where a is the orbital 
distance from the star. Once started, the evolution usually consists of three main phases. 
The first involves primarily accretion of solids onto the core, with a relatively low-mass 
envelope and a low gas accretion rate. The solids accretion rate slows down significantly 
near the point where the isolation mass (Miso) for the core is reached; for cr = 10 at 
5.2 AU this mass is about 11.5 M^. During the second phase, the gas accretion rate is 
about 3 times as high as the core accretion rate, and both are nearly constant in time. 
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Table 1: Input Parameters 






Run 


a (g/cm^) 


axY (g/cm^) Toeb (K) 


conv? 


Miso (M®) 


cjlO 


10 


700 115 


Y 


11.56 


cjlOR 


10 


700 115 


N 


11.56 


a6 


6 


420 115 


Y 


5.37 


a4 


4 


280 115 


Y 


2.92 



The envelope mass builds up relative to the core mass, which grows slowly. The third 
phase, rapid gas accretion, starts just prior to crossover mass, Afcoio = -^^onv = v2Miso- 
The calculations reported here stop at a point well beyond crossover, where dMonv/dt = 
dMiim/dt. Beyond that point, accretion up to w 1 Mjup is very fast, and the core mass 
does n ot change appreciably. Examples of this final phase are given in iLissauer et al.l 
( 20091 ): note especially Cases 21 and 21J, calculated with cr = 10 g cm~^. 

Three main cases were calculated, with a — 10, 6, and 4 g cm~^ and with all other 
parameters identical. These will be referred to as runs crlO, ct6, and (t4, respectively. A 
fourth case (crlOR) was calculated with ct = 10 g cm"^ but without convective effects 
included in the grain calculation, although they were included in the ac tual planetary 
model calculation. The most closely analogous comparison case from iLissauer et al. 
( 2009n was their C ase 21J with a = 10 g cm~^ and the same outer boundary conditions as 



used here; however ILissauer et al.l simply set the grain opacity equal to 2% of interstellar 
values. That calculation reached a crossover mass of 16.16 M® at a time of 2.3 Myr, 
and a final total mass of 1 Mjup at 2.6 Myr. The input parameters of the four runs are 
shown in Table 1, which includes a, the gas surface density axY, the surface boundary 
temperature Tnch, a column "conv?", which indicates whether convection was included 
in the grain calculation or not, and the isolation mass Miso. 

The results of our calculations with convection are shown in Figs. |T] and JH For the 
CT = 10 case, the formation of the planet is significantly speeded up in comparison to 
case 21J. The early core accretion phase is almost identical to that in case 21J as the 
envelope opacity has little effect on the rate of buildup of the core. The first luminosity 
peak is at a level of 10^^ L© at 0.4 Myr; the isolation mass is approached at about 0.45 
Myr. However Phase 2, during which core and envelope accretion rates are relatively low, 
is significantly shorter in the present calculation. The main effect of the grain settling 
calculation is to reduce the average opacity in the envelope, allowing it to lose heat more 
readily, resulting in a higher rate of envelope contraction and therefore needing a higher 
gas accretion rate to keep the outer boundary at Ra- The crossover mass is almost 
identical to that in case 21J, but it is reached at a time of only 0.97 Myr. Beyond that 
point the gas accretion rate increases rapidly, reaching the limiting rate of 2.7 x 10^^ 
Me yr "1 at 1.03 x 10^ yr with Mcorc = 16.8 M® and M^n-v = 56.8 M®. By way of 
comparison, case 21J reached the same dA/um/di with Mcorc = 17.0 M® and A/onv = 55.5 
Mq, but the time was 2.45 x 10^ yr. 

Thus the main results are tha t, as in previous work where we examine the effect 



of opacity on planetary formation ([Pollack et al.l . 119961 : iHubickvi et al.l . 120051 ). the final 



core mass is very insensitive even to large changes in the opacity. Furthermore, compared 
to calculations with full interstellar grain opacity, not including grain coagulation and 
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Figure 1: Mass of the protoplanet as a function of time for three cases that include grain settling and 
coagulation. The solid lines denote the mass of the core, the dotted lines the mass of the H/He envelope, 
and the dot-dashed lines the total mass. All cases plotted include convection in the grain calculations. 
The assumed solid surface density a is indicated for each set of curves. 
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Figure 2: The protoplanet's luminosity as a function of time for the same cases as shown in Fig. ^ 
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Figure 3; Structure of the alO case at an age of 4.9 X 10"" yr, when the core mass is 12.07 Mgj and the 
envelope mass is 1.53 M^g. Convective and radiative zones are indicated. The total radius is 2.8 X 10^^ 
cm. The temperature at the core/envelope interface is 1.33 xlO* K. 



settling ( Hubickvi et all |2005[) , the formation time for Jupiter is reduced by a factor 6. 
The main contributor to this result is the sharp reduction in the opacity that appears 
in the outer radiative zone at temperatures ~ 500-600 K. This feature is visible in the 
structure plots shown in Figs. [3] and 01 Interior to this minimum, molecular opacities 
rather than grain opacities dominate, and the opacity increases with temperature. 

The properties of Runs crlO, crlOR, cr6, and (t4 at particular times in the evolution are 
summarized in Table 2. Here Mz — Mcorc and Mxy = -^^cnv Case crlOR is not included 
in the figures, because, as shown in Table 2, the results are very similar to those of 
crlO. The small differences in the two cases are mainly a result of numerical fluctuations 
in the calculations. The properties are indicated at the following times: the maximum 
in luminosity that occurs during the main core accretion phase, the transition between 
Phase 1 and Phase 2 where the isolation mass of the core is approximately reached, a 
time in the middle of Phase 2, the time at which the crossover mass (Aicorc — ^cnv) is 
attained, and the time when the limiting gas accretion rate is reached. Following this last 
time, the planet quickly reach es its final mas s throu gh rapid gas accretion; the details of 
this final phase are treated in lLissauer et alJ (J2009I) . 



Case cr6 is calculated with cr = 6, and as expected the formation time is longer, 
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Figure 4: Structure of the crlO case at crossover, age 9.68 X 10^ yr with core mass 16 M^. Convective and 
radiative zones are indicated. The total radius is 4.19 X 10^^ cm. The temperature at the core/envelope 
interface is 1.86 xlO* K. 
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Table 2: Results 










(tIO 


crlOR 


ct6 


ct4 




Time 


0.408 


0.405 


0.654 


0.888 




Mz 


8.53 


8.12 


3.98 


2.12 


First 


MxY 


0.064 


0.044 


0.0061 


0.0014 


Luminosity 


Mz 


1.3 X 10"" 


1.25 X 10"" 


2.3 X lO"'""' 


7.1 X 10"" 


Peak 


logL 


-4.85 


-4.89 


-5.84 


-6.53 




Rp 


28.8 


27.7 


17.3 


9.86 




Time 


0.451 


0.454 


0.838 


1.19 




Mz 


11.5 


11.5 


5.60 


2.97 


End of 
Phase 1 


MxY 


0.67 


0.61 


0.63 


0.18 


Mz 
logL 


2.5 X lO"-""' 
-5.46 


2.1 X lO"'""' 
-5.52 


1.3 X 10"" 
-6.64 


4.9 X 10"^ 
-7.49 




Rp 


34.7 


35.2 


23.1 


15.7 




Time 


0.775 


0.781 


1.295 


2.738 




Mz 


14.0 


14.0 


6.50 


3.54 


Mid 
Phase 2 


MxY 


7.0 


7.0 


3.25 


1.77 


Mz 


7.9 X 10"" 


1.2 X lO"-""' 


2.1 X 10"" 


4.4 X 10"^ 


MxY 


2.7 X lO"'""' 


3.8 X lO"'""' 


7.5 X 10"" 


1.6 X 10"" 




logL 


-5.88 


-5.77 


-6.65 


-7.52 




Rp 


47.4 


48.0 


30.8 


21.8 




Time 


0.968 


0.947 


1.625 


3.618 




Mcross 


16.09 


16.11 


7.50 


4.09 


Crossover 
Point 


Mz 


2.1 X lO"-""' 


1.8 X lO"'""' 


4.4 X 10"" 


9.9 X 10"^ 


MxY 


1.8 X 10"" 


6.1 X lO"'""' 


2.2 X lO"'""' 


4.2 X 10"" 


logL 


-5.54 


-5.52 


-6.27 


-7.13 




Rp 


57.8 


60.1 


38.6 


25.5 




Time 


1.003 


0.997 


1.865 


4.038 


Onset of 


Mz 


16.8 


17.3 


8.92 


4.74 


limited 


MxY 


56.8 


52.1 


54.3 


34.0 


gas 


MxY 


2.7 X 10"" 


2.9 X 10"" 


1.84 X 10"" 


1.51 X 10"" 


accretion 













Time is in units of millions of years, Myr. 

Mass {Mz and Mxy) is in units of Earth's mass, Mg. 

The accretion rate (Mz and Mxy) is in units of Earth masses per year, Mg 

Luminosity (L) is in units of solar luminosity, Lq . 

Radius (Rp) is in units of Jupiter's present equatorial radius, Rj. 
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but still less than that for the calculation done by iHubickvi et al.l (J2005I) for this level 



of (T. The final core mass is, however, significantly reduced, to 8.92 M®. The first 
luminosity peak is reached in 0.65 Myr, somewhat later than in case ulO. The core 
accretion rate is of course slower in case cr6, but the available planetesimal supply is 
considerably less, a compensating effect. The maximum luminosity is typically reached 
when (T has been reduced (by accretion onto the planet) to 18% of its initial value. For 
case (t6, the maximum is a factor 10 less than in case crlO. Phase 2 is reached at a time 
of 8.4 X 10^ yr with a core mass of 5.6 M®, half th at of case glO at th e corresponding 
epoch. Mainly because of the reduced core mass fsee lPollack et al.l . ll996[) . the luminosity 



is considerably lower in Phase 2 than it is for case crlO, therefore energy is released less 
rapidly, contraction is slower, and the rate of gas accretion is reduced. Thus Phase 2 is 
lengthened, but its duration remains less than 1 Myr, giving a total formation time of 
1.86 Myr. 

The case cr4 represents a disk with a — A g cm^^ at 5.2 AU. This value is only slightly 
greater than that of the minimum mass solar nebula, but note that our planetesimal 
accretion prescription neglects expulsion of solids from the planet's accretion zone, so a 
somewhat higher a is probably required in order to reproduce the core accretion rate 
calculated in this case. The first luminosity peak occurs at about the same time as in cr6, 
but the value of the luminosity is a factor of 5 lower. Crossover is reached at a core mass 
of only 4 Mg and a luminosity of 10~^ L©. As shown in Figs.[5]and|ni the temperatures 
at a given density in the envelope structure are considerably lower than in case ctIO. 
The planet's internal structure, both in Phase 2 and at crossover, is characterized by an 
inner convection zone and an outer radiative zone, with the boundary at a temperature 
of only 200 K. The deep minimum in opacity just outside this boundary allows for a 
high rate of envelope radiation and the very rapid addition of mass. The onset of limited 
gas accretion is reached after a time of 4.04 Myr, approximately the mean lifetime as 
observed for protostellar disks. The heavy element core has a mass of only 4.74 M^. 

5. Discussion and Conclusions 

The formation phase of Jupiter, as determined by the concurrent accretion of solid 
particles and gas, has been calculated under the following major assumptions: (1) The 
protoplanet remains fixed in position at 5.2 AU. (2) The planet is assumed to exhibit 
runaway accretion of solids within a disk of 100 km planetesimals. The planetesimals are 
well mixed within the feeding zone of the planet; thus their surface density is constant 
in space but varies in time as the planet accretes. Migration of planetesimals into or 
out of the feeding zone is not considered. (3) The protoplanet is assumed to be the 
sole dominant mass in the region; other embryos which might compete for planetesimals 
are not considered. (4) The solids accretion rate is determined from Safronov's formula 
(e g. 1) with gravitational enhance ment factor obtained from the 3-body numerical fits 
of iGreenzweig and Lissaueij ( 19921 ). (5) All of the planetesimal material accreted by the 



protoplanet is assumed to fall into the core; thus what is referred to as 'core mass' in 
this paper is really the total heavy element excess over solar abundances. (6) The gas 
accretion rate is determined by the requirement that the outer boundary of the proto- 
planet remain at the modified ac cretion radius (eg. 2 ) . Three-dimensional simulations 
of a planet embedded in a disk ( Lissauer et al.l . l2009f) determine in part the region of 



a planet's gravitational influence. (7) The opacity caused by grains in the envelope of 
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Figure 5: Structure of the cr4 case during phase 2 at an age of 1.768 X 10^ yr, when the core mass is 
3.21 Mfg and the envelope mass is 0.74 Mq. The convective-radiative boundary is indicated. The total 
radius is 1.20 X 10^^ cm. The temperature at the core/envelope interface is 6230 K. 
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Figure 6: Structure of the crA case at crossover, age 3.618 X 10^ yr with core mass only 4.1 M^. The 
convective-radiative boundary is indicated. The total radius is 1.93 X 10^^ cm. The temperature at the 
core/envelope interface is 7680 K. 
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the protoplanet is determined by a detailed calculation of grain settling and coagulation, 
with the grain size distribution and the resulting grain opacity being recalculated at every 
level in the envelope at every time step. It is this last assumption that represents a major 
improvement over past calculations of Jupiter's formation, and it results in a significant 
reduction in the formation time scale as compared to past calculations. There are ap- 
proximations involved in the treatment of convection zones in the grain calculations, but 
the comparison of cases crlO and crlOR shows that these have a negligible effect on the 
results. The grains in the incoming gas are composed of silicates only, with an abundanc e 
by mass of 1% that of the gas. Molecular opacities are included ( Freedman et al.l |2008[) . 



The main parameter varied in this set of calculations is the solid surface density a 
at 5.2 AU in the primordial disk; values considered are 10, 6, an d 4 g cm~^. A s in 
previous simulations of Jupiter's formation using the prescription of lLissaueii ( 19871 ) for 



the runaway accretion of solids, the results show that there are three main phases. Phase 
1 involves primarily the accretion of solid particles to build up the core; the mass of the 
gaseous envelope is fairly negligible. Phase 2 begins when the isolation mass is reached; 
accretion of gas takes place at about three times the rate of accretion of solid particles, 
and the envelope builds up in mass relative to the core. Phase 3, rapid gas accretion, 
begins close to the point where core and envelope masses are equal. The calculations in 
this paper are taken beyond the crossover mass up to the point where the gas accretion 
rate becomes limited by the rate at which the disk can supply the mass, as determined 
by three-dimensional simulations (1.5 — 3 x 10^^ M^ yr~^). The accretion up to the final 
mass of Jupiter takes place quickly after that point ( Lissauer et al.l . 120091 ) . 



The formation time to crossover is dominated by Phase 2 in all simu lations. Case 
(tIO may be compared with a very similar calculation ( Lissauer et al.l . 120091) in which the 



grain opacity was simply assumed to be 2% that of interstellar grains. Phase 1 in the 
present case was completed in about the same time as in the comparison case, but Phase 
2 (0.52 Myr) was a factor of 3.6 shorter, giving a total formation time of about 1 Myr, 
compared with 2.6 Myr in the comparison case. Correspondingly, the mean luminosity 
in phase 2 was a factor of 3.6 higher in the present run than in the comparison case. The 
main cause of the increase in luminosity is the very low opacity in the region around 600 
K in the present model (Fig. [3]), with a minimum of opacity of about 6 x 10"'* cm^ g~* 
as compared with 0.04 cm^ g^* at the same temperature in the comparison case. Some 
layers in the present case have opacities higher than 2% interstellar, but they are mainly 
in the outermost layers that are either convective or rather optically thin, so they do 
little to influence the flow of radiation. In spite of the difference in formation time for 
the two cases, the core mass at onset of limiting gas accretion was practically the same in 
both case s, about 17 Mm. This ma ss is higher than that usually derived for the present 
Jupiter (jSaumon and Guilloti . 120041) b ut it is close to t hat pr edicted from an independent 



calculation of the equation of state bv lMilitzer et al.l (|2008l ). 



The question then arises whether a reduction in ct, which is known to result in a 
smaller core mass, still leads to a reasonable formation time. Cases cr6 and (t4 show 
that this is indeed the case. The formation times are 1.9 and 4.0 Myr, respectively. The 
ratios of Phase 2 times for the three cases scale as M.J^ /L, as derived bv IPollack et al.l 
(|l996. ) under the assumption that planetesimal accretion provides most of the luminosity 
during Phase 2. The value 4 Myr is w ithin the range of e xpected lifetimes of observed 



circumstellar disks around young stars ([Haisch et al.l . 120011 ). Note that the cr6 case has a 
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format ion time only 14% as lo ng (and Phase 2 time only 8% as long) as the corresponding 



case m 



Hubickvi et al.l (|2005l ) with grain opacities set at 2% of interstellar values, while 
the crlO case took 40% as long as the corresponding case (lOLcx)) in that earlier work. 

The main conclusions of this work are: (1) Within the basic framework of a core- 
nucleated accretion model for Jupiter's formation, the effects of grain settling and coag- 
ulation in the gaseous envelope of a protoplanet are significant in determining the gas 
accretion rate and formation time. (2) The formation time for a Jupiter at 5.2 AU with 
a solid surface density of 10 g cm~^ (about 3 times that in the minimum mass solar 
nebula) is about 1 Myr, significantly shorter than that in comparable comparison cases. 
If a is taken to be 4 g cm~^, the formation time is still reasonably short, about 4 Myr. 
(3) In the case with tr = 4 g cm~^, Jupiter's core mass i s calculated to be 4. 7 4 M 



consistent with the Jupiter structure models computed bv lSaumon and Guillod ( 200 



No arbitrary cutoff in the core accretion rate was needed to obtain this result. 

Further work is needed, however, to confirm this result. The main points which 
need to be improved include, first, the planetesimal accretion model, which should be 
replaced by a statistical calculation taking into account the evolution of the planetesimal 
size distribu tion and the effects of the p rotoplanet itself on the planetesimal velocity 



distribution (jWeidenschilling et al.L Il997f) . Numerous estimates of the Phase 1 accretion 



time have been p ublished. One example is the (statistical -I- N-body) calculation of 



Inaba et al.l ( 20031 ) which takes into account planetesimal fragmentation as well as the 



effect of the gaseous envelope in increasing the effective cross-section. A planetesimal size 
of 10 km is assumed. If the opacity in the planetary envelope is reduced by a factor 100 
compared with interstellar values, a core approaching 10 Mq is reached in 5 Myr at 5.2 
AU; however a is required to be about 5 times that in the minimum mass solar nebula. 
This and other simulations, along with the present results, suggest that the formation 
time for Jupiter may be controlled by Phase 1, not Phase 2. Second, the effects of the 
ices should be included in the calculations of grain settling and coagulation. Third, the 
effects of the dissolved accreted solid material, primarily ices, should be included in the 
equation of state and opacity of the envelope and in the calculation of the solid core 
mass. 
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